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Intermittency in fluid turbulence can be evidentiated through the analysis of Prob- 
ability Distribution Functions (PDF) for velocity fluctuations, which display a 
strong non— gaussian behavior at small scales. In this paper we investigate the oc- 
currence of intermittency in plasma turbulence by studying the departure from the 
gaussian distribution of PDF for both velocity and magnetic fluctuations. We use 
data coming from two different experiments, namely in situ satellite observations 
of the inner solar wind and turbulent fluctuations in a magnetically confined fusion 
plasma. Moreover we investigate also time intermittency observed in a simplified 
shell model which mimics 3D MHD equations. We found that the departure from a 
gaussian distribution is the main characteristic of all cases. The scaling behaviour 
of PDFs are then investigated by using two different models built up in the past 
years, in order to capture the essence of intermittency in turbulence. 



1 Introduction 

The statistics of turbulent fluid flows can be characterised by the Probability 
Distribution Function JPDF) of velocity differences 6u r — u(x + r) — u(x) over 
varying scales r (see El and references therein). At large scales the PDF is 
approximately gaussian, as the scale decreases, the wings of the distribution 
become increasingly stretched, so that large deviations from the average value 
are present. This phenomenon is usually ascribed to intermittency. The way 
the PDF departs from the gaussian is very interesting because this has to do 
with models for intermittency in turbulence. If we introduce a scaling for 
fluctuations, namely 5u r ~ r h as the fluid equations indicate, a change of scale 
r — > It (with / > 0) leads to 

Suir = I Su r 

This is interpreted as an equality in lawQ, that is the right-hand-side has the 
same statistical properties as the left -hand-side. It can be easily shown that if 
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h is unique, that is in a pure self-similar (fractal) situation, the standardized 
variables Su r / (Su 2 .) 1 ' 2 have the same PDF for all scales r. On the contrary 
me must invoke the multifractal theory (or some other models) to describe 
intermittency, by conjecturing the presence of an entire range of values of h. 

Intermittency has been observed some time ago and deeply investigated in 
fluid flowsEJ, .and more recently also in Magnetohydrodynamic (MHD) flows (see 
for examplcH). In particular, experimental studies carried out in MHD flows 
inly by analyses of satellite measurements of solar windJluctuations 
• by using high resolution 2D numerical simulations Qii3 and Shell 
Models EUO. All these analysis deal with the scaling exponents of structure 
functions, aimed to show that they follow anomalous scaling laws which can 
be compared with the usual energy cascade models for turbulence. 

The non gaussian nature aLPDF in MHD solar wind turbulence has been 
evidenced by Marsch_and Tuliia, followed by quantitative analysis made by 
Sorriso-Valvo et al. Ed. In fact, in order to investigate the properties of in- 
termittency through the analysis of non gaussian character of PDF, it would 
be necessary to quantify the departure of PDF from gaussian statistics and to 
analyse how this departure depends on the scale r. In this work we present 
two models, introduced earlier to describe intermittency in fluid flows. Both 
models describe the evolution of PDF with r, by allowing some free parameters 
to evolve with the scale r. We then analyse some data of plasma turbulence 
coming from very different situations. First of all we show that both mod- 
els are able to capture time intermittency observed in a discrete model which 
describes the MHD turbulent cascade (Section 3), then we show that also in- 
termittency in real MHD turbulence, say the turbulence observed in the Solar 
Wind (Section 4) and magnetic turbulence as observed in a fusion device (Sec- 
tion 5), is described by the models. We will describe intermittency in plasma 
turbulence through the scaling laws of the characteristic parameters of both 
models. Finally (Section 6) we discuss the results we have obtained. 



2 Models for PDF 

As stated above, in this section we briefly introduce two models which describe 
the scaling evolution of PDF, trying to describe intermittency through the scale 
variation of characteristic parameters. We consider thfi-jpodel based on a log- 
normal energy cascade proposed by Castaing et al. llj'E3; the second model 
describes the PDF through a sequence of stretched exponentials and can be 
formally derived through the theory of extreme deviationstLH. 
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2.1 The energy cascade model 

Because of the idea of self-similarity underlying the energy cascade process 
in turbulence, Castaing and co-workers LLj introduced a model which tries to 
characterize the behaviour of the PDFs through the scaling law of a parameter 
describing how the shape of the PDF changes in going towards small scaleslill 
In its simplest form the model can be introduced by saying that the PDF of the 
increments dtp (representing here both velocity and magnetic fluctuations) at 
a given scale r, is made by a convolution of the typical Gaussian distribution 
Pq{5iP,(t) — {\f 7 hto)~ 1 exp (— Si/j 2 /2a 2 ), with a function G T (er, A) which rep- 
resents the weight of the gaussian distribution characterized by the standard 
deviation a: 

P T (8i/>, A) = J G T (a, A) P G (Si/J, a) da , (1) 

where A is a parameter representing the width of G T . 

In the usual approach where the energy cascade is introduced through a 
fragmentation process, a is directly related to the local energy transfer rate e. 
In a self-similar situation, where the energy cascade generates only a scaling 
variation of a = (Sip 2 ) 1 / 2 according to the classical Kolmogorov's pictured, 
G T (o~) reduces to a Dirac function G T (a) = S(a — cfq). In this case from eq. 
(|l|) a Gaussian distribution P T (5ip) — Pg(<5^>, oo) is retained. On the contrary 
when the cascade is not strictly self-similar, the width of the distribution G T 
is different from zero. In this way the scaling behaviour of the parameter 
A (which takes into account the height of the PDFs wings) can be used to 
characterize intermittency. 

2.2 A model coming from Extreme Deviations Theory 

Recently Frisch and Sornette0, have built up a more general fragmentation 
model which is based on the extreme deviations theory (EDT), which yield 
non-gaussian PDF as the number of fragments increases. Let us consider a se- 
ries of random events Xi , independent and distributed according to a probabil- 
ity density p(x). According to the central limit theorem, the sum y = X)T=i x * 
is distributed according to a gaussian function, that is fluctuations y/n are 
0(n -1 / 2 ). On the contrary EDT deals with events where fluctuations are 
O(l), with a finite value n of fragments and large values of y. The theory al- 
lows us to recover the asymptotic shape of the PDF P n (y) starting from p(x). 
The authors 113 consider the general case where p(x) — exp[—f(x)], and f(x) 
is differentiable and increases quite rapidly as x — ► oo. Assuming that the 
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second derivative is positive / (ar) > and that x 2 f"(x) — ► oo (for large x), 
the asymptotic form for the PDF of y is given by 



Pn(y) ~ exp[-nf(x/n)] . 



The probability conditioned on Xi whose sum is y, is localized for large y in 
a "democratic" manner near the value for which y/n ~ X4 ~ a; 2 ~ ... — x n . 
If during the fragmentation process the fragments have probability p(x) ~ 
exp(— Cx 7 ) with 7 > 0, and y = logS 1 , where S 1 = w 1 w 2 ...w n (that is Xj = 
logWi), then 



From the turbulence point of view, if at the scale r the value of each 5ij) is due 
to the product of fragments derived from structures at larger scale, from EDT 
we must expect a PDF as a stretched exponential 



where /3(r) and h(t) are the parameter we must fix from data. 

With these models in mind we will discuss now what kind of results are 
obtained for plasma turbulence. 

3 A shell model for MHD turbulence 

Turbulence is a phenomenon in which chaotic dynamics and power law statis- 
tics coexist. In general, due to the high Reynolds numbers involved in realistic 
situations, direct numerical simulations of 3D turbulence is not possible. In 
order to avoid this problem and in order to capture the gross features of turbu- 
lence, an alternative approach can he followed by investigating reduced models 
of turbulence, such as shell modelsLL3. In these models the wave vector domain 
is spaced exponentially through k n = fco2" (n — 1,...,N) and the nonlin- 
ear dynamics is reproduced by using a single complex variable for each n-th 
shell, that is u n (t) for the velocity field and b n (t) for the magnetic field. The 
quadratic nonlinear interactions between modes are then recovered by impos- 
ing that: i) the dynamical equations for each shell variable must contain the 
same quadratic terms of the original equations; ii) the quadratic invariants of 
equations must be correctly preserved in the inviscid case; iii) the nonlinear 
coupling occur between nearest or next nearest shells, modeling local interac- 
tions in the k space. By using these criteria an MHD shell model whicruis the 
.analogous of the GOY model for hydrodynamics is easily obtained (see t3 and 
113 for a complete review on shell models): 



P(S) 



exp(-CnS' 7 /") . 



P(6i/)) ~exp -/3(r)|(5Vr (T) 



(2) 
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dUn 

dt 



(u n+ iU n+2 - b n+1 b n+2 ) - -(tt n -l«n+l - b n -ib n+ i) + (3) 



-(u„_ 2 Un-l - &n-2&n-l) 



- f^W n + /„ 
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7-(«n+l&n+2 - 6 n+l fi n+2) + \{u n -\b n+ i ~ b n -iU n+1 )+ (4) 

6 6 



+ ^(«n-2&n-l - &n-2ttn-l) 




where ^ and r? are respectively the kinematic viscosity and the resistivity, and 
f n and g n are forcing terms. Equations (^|) and (|^), in the inviscid and unforced 
case, conserve the total energy 



the cross-helicity 



1 N 



i N 



and the magnetic helicity 



N 



ff(i) = £("l 



Apart for other features, shell models share with original equations scaling 
laws of the structure functions S p (n) = (\x n \ p ) ~ k n where x n is either u n 
and b n . Scaling exponents Cp deviates from the linear Kolmogorov scaling as 
a consequence of the intermittency in the dynamical system. 

We solved numerically the system of equations for the shell model, by using 
the following set of parameters: N = 23 shells, v = rj = 10 -9 , the integration 
step dt — 5 x 10~ 5 . The forcing term acts only on the velocity variables (i.e. 
g n = 0) at large scales (n = 1,2), the magnetic field being generated by a kind 
of dynamo effect which is at work in the MHD shell model £3. The forcing, 
which is assumed to be random, is calculated according to a Langevin equation 



5 



df n /dt = —f n TQ+9, where tq is the characteristic time of the largest scales, and 
9 is a gaussian white-noise. PDFs for both the velocity and magnetic variables 
have been obtained and are shown in figures [l] and || for different shells. As it 
can be seen PDFs depends on the scale l n ~ k~ x . They are almost Gaussian 
at the large scales (small shells k n ) and are increasingly non-Gaussian at small 
scales (large shells k n ). 

In order to make a quantitative analysis of the energy cascade leading to 
the process just described, we have tried to fit the distributions by using, for 
the model ([[]), the log-normal ansatzEj 

lnV/ao) 



G T (a) da = == exp 

V ' A(t)v^F 



2A 2 (t) 



dQna) (5) 



even if also other functions gives rise to results not really different. The pa- 
rameter (To represents the most probable value of a, while A 2 (r) = ((Alncr) 2 ) 
is the width of the log-normal distribution of a. 

We have fitted the expression (0) on the PDFs obtained in the shell model 
for both velocity and magnetic intensity, and we have obtained the correspond- 
ing values of the parameter A 2 . The values of the parameters <7o, which do not 
display almost any variation with £ n are reported in table I. In figure |l| we 
plot, as full lines, the curves relative to the fit. As can be seen the scaling 
behaviour of PDFs in all cases is very well described by the function (Q). From 
the fit, at each scale £ n , we get a value for the parameter A 2 (£ n ), and in figure 
^ we report the scaling behaviour of A 2 (£„). 

The scaling evolution of PDFs have been analysed also using the stretched 
exponential model P(8xi) ~ exp[— /3(&c^) M ]; the results of the fit are displayed 
in figure ^| as a full line. The parameters fi found from the fit are shown in 
figure [| as a function of the scale £„, for both xi n — u n and xg n = b n . As 
can be seen, a power law [i(£ n ) ~ £® is found in a very wide range for the 
magnetic field; the range is less wide in the case of the velocity field, for wich 
a saturation of fj, is reached. The minima values for /i are reported in the table 
||, together with the exponents a found for the power laws. 

As we will show in the next two Sections, both models are able to capture 
intcrmittency also for turbulence in non homogeneous situations, namely the 
Solar Wind and a fusion plasma. To avoid superpositions we describe for 
the Solar Wind only the results obtained with log-normal ansatz, while for 
the fusion plasma we describe only the results obtained with the stretched 
exponentials, even if in both [l] and ^| the reader can find the characteristic 
parameters obtained through the fit of both models for both samples. 
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Figure 1: The PDFs of the coefficients b n and u n for three different values of n. The fit with 
the log-normal model by Castaing et al. QJLT) is superimposed as a full line. 
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Figure 2: The values of the parameter A 2 as a function of the scale l/ta. Black circles are 
for the magnetic field, while the white circles are for the velocity. 
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Figure 3: The PDFs of the coefficients b n and u n for three different values of n. The fit with 
the stretched exponentials functions is superimposed as a full line. 




Figure 4: The values of the exponent /i as a function of the scale £/£o- Black circles are for 
the magnetic field, while the white circles are for the velocity. 
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4 Solar Wind Observations 

The satellite observations of both velocity and magnetic field in the interplan- 
etary space, offer us an almost unique possibility to gain information on the 
turbulent MHD state in a very large scale range, say from 1 AU (Astronomical 
Unit) down to 10 3 km. Here we limit to analyse only plasma measurements of 
the bulk velocity V(t) and magnetic field intensity B(t). We based our analy- 
sis on plasma measurements as recorded by the instruments on board Helios 2 
during its primary mission in the inner heliosphere. The analysis period refers 
to the first four months of 1976 when the spacecraft orbited from 1 AU, on 
day 17, to 0.29 AU on day 108. The original data were collected in 81s bins 
and we choose a set of subintervals of two days each. The subintervals were 
selected separately within low speed regions and high-speed regions selected 
in a standard way according to a threshold velocity til. For each subinterval 
we calculated the velocity and magnetic increments at a given scale t through 
SV T = V{t + r) — V{t) and SB T = B(t + r) — B(i), which represent characteris- 
tic fluctuations across eddies at the scale r. Then we normalize each variable 
to the standard deviation within each subinterval Sv T = SV T /(SV 2 ) 1 ^ 2 and 
5b T = 8B T / {5B 2 .) 1 / 2 (where brackets represent average within each subinterval 
at the scale r). We calculate the PDFs at ff different scales logarithmically 
spaced r = At 2™, where n = 0, 1, 10 and At = 81s. We collect the num- 
ber of events within each bins by using 31 bins equally spaced in the range 
within 3 times the standard deviation of the total sample. Before we mixed 
the different subperiods belonging to a given class (high or low speed streams), 
we tested for the fact that the gross features of PDFs shape does not change 
in different subintervals. Then our results for high and low speed streams are 
representative of what happens at the PDFs. 

The results are shown in figures |s| and ||, where we report the PDFs of 
both velocity and magnetic intensity for both the high speed streams and 
the slow speed streams. At large scales the PDFs are almost Gaussian, and 
the wings of the distributions grow up as the scale becomes smaller. This is 
true in all cases, say for both types of wind. Stronger events at small scales 
have a probability of occurrence greater than that they would have if they 
were distributed according to a gaussian function. This behaviour is common 
for intprmittpnry as currently observed in fluid flows u and in the solar wind 
turbulence B0. 

A quantitative_analysis of the energy cascade has been made by using the 
log-normal ansatzlij rta) . We have fitted the expression (j|) on the experimental 
PDFs for both velocity and magnetic intensity, and we have obtained the 
corresponding values of the parameter A 2 . The values of the parameters <jq, 
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Figure 5: The scaling behaviour of the PDF for 5v T as calculated from the experimental 
data (white symbols) in the fast (left) and slow (right) streams. The full lines represent the 
fit obtained through the model by Castaing et al. (Ill) . as described in the text. 
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Figure 6: The scaling behaviour of the PDF for 8b T as calculated from the experimental data 
(white symbols) in the fast streams. Again, the full lines represent the fit obtained through 
the model by Castaing et al. (hi), as described in the text. 
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Figure 7: We show the scaling behavior of A 2 (r) vs. t for both fast (black symbols) and 
slow (open symbols) streams. Circles refer to the magnetic field intensity, squares refer to 
the bulk velocity. 



which again do not display almost any variation with r are reported in table 
1. In figures || and |(] we plot, as full lines, the curves relative to the fit. As can 
be seen the scaling behaviour of PDFs in all cases is very well described by the 
function (Q). From the fit, at each scale r, we get a value for the parameter 
A 2 (t) , and in figure we report the scaling behaviour of A 2 (r) for both high 
and low speed streams. Starting from a very low value at the scales of about 1 
day, A 2 increases abruptly to A 2 ~ 10 _1 at about 2 hours, and finally a scaling 
law starts to become evident up to At = 81 sec. In this last range, which 
corresponds roughly to what is usually called the "Alfvenic range" , we fitted 
the parameter with a power law A 2 (r) ~ r _Q . The values of a obtained in 
the fitting procedure and the corresponding range of scales, are reported in 
the table |l|. In the same table, we report the values of the maximum value of 
A 2 (t) in the range (Afrti, namely A 2 ^^., which is an indication for the strenght 
of the intemittencyliSEi. 

Looking at figure and at A 2 raaa . , it can be seen that both in fast and in 
slow streams magnetic field intensity is more intermittent than bulk velocity 
(values of A 2 are at least two times larger for magnetic field intensity than 
foil velocity). The same indications comes from 2D MHD xLirect simulations 
□113, and in different analysis of solar wind intermittency £3. The values of 
A 2 (t) are more or less the same for magnetic field intensity both in fast and 
in slow wind. This is due to the fact that magnetic field intensity fluctuations 
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are related to compressive fluctuations, which should have the same nature in 
both types of wind. The bulk velocity fluctuations on the contrary are more 
intermittent at small scales (81 sec) in the fast wind and at large scale (~ 1 
hour) in the slow wind. From table [l] it appears that the palue of a is not 
universal, a result which has also been found in fluid flows E3. The fact that 
the value of a for the velocity fluctuations in slow wind is the same as the 
value of a for magnetic field intensity suggests that intermittent structures in 
slow wind are perhaps mainly associated with compressive fluctuations. On 
the contrary the different value of a found in fast wind evidences a different 
nature of velocity fluctuations in fast wind, perhaps related to the fact that 
such fluctuations are mainly incompressible. This result is in agreement with 
what has been found by Veltri and MangeneyE3. These authors found that in 
fast wind the more intermittent structures are tangential discontinuities with 
almost no variation in magnetic field intensity, while in slow wind the most 
intermittent structures are shock waves, which display the same behaviour in 
bulk velocity and magnetic field intensity. 

5 The Reverse Field Pinch: RFX experiment 

In this section we report evidences for the presence of intermittency in another 
type of magnetized fluid, namely a plasma of interest for controlled thermonu- 
clear fusion research, confined in reversed field pinch (RFP) configuration. 

The RFP is a toroidal configuration of magnetic fields which is pxoposed 
as an alternative to the tokamak for confining fusion grade plasmas E3. The 
configuration is characterized by toroidal and poloidal magnetic field compo- 
nents of comparable magnitude (in a tokamak the field is mainly toroidal). 
The configuration is a near-minimum energy state to which a plasma relaxes 
under proper constraints Ej. The toroidal field changes sign in the outer part 
of the plasma, a feature which gives the name to the configuration. Such field 
reversal, which improves the MHD stability of the configuration, is sponta- 
neously generated by the plasma, and is maintained against resistive diffusion 
by the dynamo process □. This is achieved through the action and nonlinear 
coupling of several resistive magnetohydrodynamic modes, which give rise to 
a high level of magnetic turbulence (of the order of 1% of the average field 
in present day experiments, i.e. two orders of magnitude larger than in toka- 
maks). This high fluctuation level makes the RFP very suited for the study 
of MHD turbulence properties, mainly for their magnetic part. The magnetic 
turbulence has been demonstrated to be the main cause of energy and particle 
transport in the RFP core, whereas at the edge its contribution is still under 
investigation. In this region the electrostatic turbulence has been proved to 
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Figure 8: The scaling behaviour of the PDF of the magnetic fluctuation Sb T measured in the 
RFX experiment. The full lines represent the fit obtained through the stretched exponential 
model, as described in the text. 



give an important contribution to the particle transport El 

The measurements have been obtained in the RFX experiment, which is 
thc-Largest RFP presently in operation (major radius 2 m, minor radius 0.457 
m) c3. RFX is designed to reach a plasma current of 2 MA, and currents up 
to 1 MA have been obtained up to now. The measurements were performed 
in low currents discharges (300 kA) using a magnetic probe inserted in the 
edge plasma. The probe consists of a coil housed in a boron nitride measuring 
head. The coil measures the time derivative dtB of the radial component B(t) 
of the magnetic field. The radial direction in this case goes from the core 
plasma to the edge. The sampling frequency of the measurements is 2 MHz. 
Measurements have been collected at different values of the probe insertion 
X = r/a into the plasma, X — being the position of the inner convolution of 
the RFX graphite first wall, r being the probe radial position and a being the 
radius of the ring. Different measurements arc also identified through a shot 
number. In RFX two different components of the magnetic fluctuations can 
be identified: a localised and stationary magnetic perturbation, originated by 
the tearing modes responsible for the dynamo which tend to be phase-locked 
and locked to the wallEZI, and an high frequency broadband activity, which 
is investigated here. All measurements presented were made away from the 
stationary perturbation. 
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Figure 9: We show the scaling of c(r) vs. r for the RFX magnetic field at different insertion 
values: X = 0.97 (black circles), X = 0.95 (white circles) and X = 0.91 (stars). 



We calculated the characteristic magnetic fluctuations at the scale r, namely 
8B(t) — B(t + t) — B(t). For each position within the device, we can study 
the statistical behaviour of fluctuations at different scales r. In figure || we 
report the PDFs of Sb T = SB r / (SB 2 ) 1 / 2 at different scales for a given value of 
X. As it can be seen the PDFs do not collapse to a single curve, but follow 
the usual characteristic scaling behaviour evidenced also in solar wind mea- 
surements. This behaviour is visible for all the values of the insertion X. A 
fit to the PDFs we obtained has been made by using a stretched exponential, 
namely P(Sb T ) ~ exp[— /3((56 T ) /i ], and their scaling behaviour is then inves- 
tigated by looking at the scaling laws of both the parameters (3(t,X) and 
/i(r, X). The parameter fi follows a power law in a certain range of scales, that 
is: (i(r,X) ~ T a ( x \ even near the external wall (figure 5). 

The scaling exponent a(X) (see table ||) remains near to the value 0.1 in 
the different shots, and the range of scales where a scaling for the parameter 
/z(t, X) is visible is larger for measurements near the first wall with respect to 
measurements inside the device. The "strenght" of the intermittency effects 
can be compared for the different insertion values X, by comparing the min- 
imum values found for the exponent fJ,(r) in each shot. As can be seen from 
the table (||), such values increase going toward the core, showing that the 
intermittency is stronger near the wallet. 
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Table 1: We report the values of the parameters <ro, ^max an d a obtained in the fitting 
procedure for A 2 (t) of the log-normal model in the shell model data, solar wind and RFX 
(shot 8414). We also report the range of scales where the fit has been done. For the shell 
model, the range of fit is in I/Iq units. 

g o tfnax a Scales (SW: hours) 

u SM 0.83 ±0.04 0.44 ±0.10 0.15 ±0.07 1.5 x 10" 5 + 10~ 3 

b SM 0.79 ±0.04 0.40 ±0.10 0.19 ± 0.06 3.8 x 10~ 6 + 10~ 3 

Bsw(fast) 0.85 ± 0.05 0.88 ±0.04 0.19 ± 0.02 0.02-^0.72 

Bsw(slow) 0.90 ± 0.05 0.73 ±0.04 0.18 ±0.03 0.02-^0.72 

Vswifast) 0.90 ± 0.05 0.51 ± 0.04 0.44 ±0.05 0.02 -e- 1.44 

Vsw(slow) 0.95 ± 0.05 0.37 ±0.03 0.20 ±0.04 0.02 -e- 1.44 

RFX siu 1.01 ± 0.02 0.21 ±0.01 0.42 ±0.03 l^s^-36^s 



Table 2: The values of Umin and the scaling exponents a, obtained for the scaling of the 
exponents fj, of the stretched exponential fit of the data, from shell model, the RFX and the 
solar wind. The scale ranges in which a power law has been found are also shown. In the 
RFX, different shots atdiffcrent insertions X = r/a are displayed. Two inner shots (8425 
and 8426) do not display clear scaling laws. However, we report the minima of fi. For the 
shell model, the range of fit is in £/£q units. 





X 


f-^min 




a 


Scales (RFX: /us) 


USM 




0.97 ±0.03 


0. 


.04 ±0.01 


1.5 x 10~ 5 -^3.9x 10 


bsM 




0.99 ±0.03 


0. 


.04 ±0.01 


3.8 x 10~ 6 -^6.3 x 10 


RFXs414 


0.97 


1.19 ±0.02 


0. 


.10 ±0.01 


l-e-36 


RFX S 417 


0.96 


1.25 ±0.01 


0. 


.09 ±0.01 


l-e-128 


RFX$420 


0.91 


1.23 ±0.01 


0. 


.10 ±0.01 


8-^ 128 


RFX$425 


0.90 


1.50 ±0.02 








RFXg426 


0.86 


1.55 ±0.03 








Bsw(fast) 




0.74 ±0.01 


0. 


.16 ±0.01 


0.02ft 24ft 


Bsw{slow) 




0.81 ±0.02 


0. 


.16 ±0.02 


0.02ft + 24ft 


Vswifast) 




0.92 ±0.02 


0. 


.12 ±0.02 


0.02ft + 24ft 


Vsw(slow) 




1.04 ±0.03 


0. 


11 ±0.01 


0.02ft + 24ft 
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6 Discussion and conclusions 



The departure from gaussianity of the PDFs of fields increments is well known 
in fluid turbulence, and has been pointed out in MHD plasmas. The data 
presented in this paper clearly confirm this peculiar behavior, showing that 
intermittency strongly affects the scaling behavior of MHD turbulence. The 
analysis of data with the two models presented here has a duplex interest: i) 
it is an experimental test for the theoretical models; ii) it gives a quantitative 
characterization of intermittency in different physical systems, leading to the 
the possibility of a comparison between them. 

For the first point, we showed and discussed in detail here the results 
obtained when we used the log-normal model to describe intermittency in 
the solar wind turbulence, and the stretched exponential model to describe 
intermittency in RFX plasma. However it can be shown that both models are 
able to capture the scaling behavior of PDFs in both the solar wind and RFX 
experiments, just as in the case of the MHD shell model. Very gpod fits are 
obtained, and very close values of the x 2 are found in both cases c3. In table 
|l| we report the results relative to the fit of PDFs obtained in RFX plasma 
(namely the shot number 8414) with the model by Castaing et al. Ii3. Other 
samples give similar results. In table ^| we report the results obtained when 
we fit the PDFs calculated with the dataset of the solar wind turbulence, with 
the stretched exponential model. 

Concerning the second point, we first observe that the two models give con- 
sistent results. The similar role played by X 2 max m the log-normal model and 
fJ-min m the stretched exponential model, as well as the slopes of the power-low 
scaling of A 2 and fj,, is evident. We can thus look at such parameters to com- 
pare the intermittency of the magnetic field in both experimental datasets we 
analysed. As a first comparison we can observe, looking directely at the PDFs 
shape of magnetic fluctuations, that intermittency in solar wind is stronger 
than in the RFX experiment plasma. The values of both parameters X^ax an d 
Mmm ( see tables [j] and ^) confirm such result. Clear scaling laws are found for 
the parameters of both models, so that it seems difficult to prefer one model 
with respect to the other on the basis of our analysis. 

Even if the experimental results we presented do not give indications 
wheter one model is better than the other in fitting experimental results, the 
stretched exponential is a more flexible function with respect to (Q). This 
is only a "technical" consideration and does not involve a theoretical prefer- 
ence. From our experience the code developed for the fit with the stretched 
exponential is more performant then the code developed by using the log- 
normal model. It is however worthwhile to compare intermittency found in 



1G 



real datasets with time intermittency found in the simplified models for MHD 
turbulence. Looking at the PDFs in figures [l], [| and we can realise that 
intermittency in the shell model is similar to that obtained in the solar wind, 
that is turbulence in the shell model is more intermittent than turbulence in 
RFX plasma. Furthermore looking at the values of fA m in, it can be seen that 
magnetic turbulence in the solar wind is more intermittent than in the shell 
model, while for the velocity the intermittency is comparable in both cases. 

To conclude we would like Id remark that numerical simulations of MHD 
equations in 2D configurations E£l display strong intermittency which has been 
successfully described by the model introduced by Castaing et al. t3. 
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